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; ABSTRACT 

We study the possibility of detecting oscillating patterns in the equation of state 

I (EoS) of the dark energy using different cosmological datasets. We follow a phenomeno- 

logical approach and study three different oscillating models for the EoS, one of them 
periodic and the other two damped (proposed here for the first time). All the models 
are characterised by the amplitude value, the centre and the frequency of oscillations. 
In contrast to previous works in the literature, we do not fix the value of the frequency 

\ to a fiducial value related to the time extension of chosen datasets, but consider a dis- 

crete set of values, so to avoid arbitrariness and try and detect any possible time period 
in the EoS. We test the models using a recent collection of SNcIa, direct Hubble data 

\ and Gamma Ray Bursts data. Main results are: I. even if constraints on the amplitude 

are not too strong, we detect a trend of it versus the frequency, i.e. decreasing (and 
even negatives) amplitudes for higher frequencies; II. the centre of oscillation (which 
corresponds to the present value of the EoS parameter) is very well constrained, phan- 
tom behaviour is excluded at Itr level and trend which is in agreement with the one 
for the amplitude appears; III. the frequency is hard to constrain, showing similar sta- 
tistical validity for all the values of the discrete set chosen, but the best fit of all the 
scenarios considered is associated with a period which is in the redshift range depicted 

\ by our cosmological data. The "best" oscillating models are compared with ACDM 

using dimensionally consistent a Bayesian approach based information criterion and 
the conclusion reached is the non existence of significant evidence against dark energy 

\ oscillations. 

in the framework of a self-consistent theoretical cosmological 
model? This is the main task of modern cosmology and no 
unique answers have been given so far. 

The ACDM model is the simplest (from a statistical 
point of view) and the most accepted (it is called concor- 
dance model just for this reason), and in this scenario ac- 
celeration is driven by the famous cosmological constant. A, 
which contributes to the energy /matter content of the Uni- 
verse by more than a 74%. At the same time it requires 
the presence of a large amount of cold dark matter (about 
the 22% of total energy/matter content), i.e. non-baryonic 
matter which does not interact with electromagnetic radia- 
tion, but it is detectable only by its gravitational interaction 
with ordinary baryonic matter (which makes the remaining 
4%). This mode l prov i des a good fi t to m o st of the data 
( Tegmark erakl (|2004l ): ISeliak et all (|2005l ): ISanchez et al] 
giving a reliable snapshot of the current Universe; 
but it is also affected by serious theoretical shortcomings 
that have motivated the search for alternative candidates 
generically referred to as dark energy or quintessence. Such 
models range from conventional scalar fields rolling down 
self-interaction potentials, to non-canonical scalar field mod- 
els (phantom, k-essence,etc.); from phenomenological unified 



1 INTRODUCTION 



' We have nowadays a great amount of independent data 
' sets available for studying the present dynamical state of 
the Universe. High quality data com ing from th e Hub - 
ble diagram of Ty pe la Supernovae jRiess et al.l (|2004h : 
lAstier et all (|200d ): IClocchiati et all (|2006l ')'): the mea- 



surements of cluster properties as the mass, the cor- 
relation function and the evolut io n with red s hift o f 
their abundan ce fEkc et al.' ^199^); IViana et all (|2002l ): 



iBahcall et aTl ((2003'): Bali caU fc Bode] 1200^1 1: t h e opt 
cal surveys of l a rge scale structure |Pope et "ahl (|2005h ; 
ICole et all (|2005l '): lEisenstein et all (|2005l Vh the anisotropies 
of the cosmic microwave b ackground (|de Bernardis et al.l 
(|2000h : ISpergel et all (|2003l ')'): the cosmic shear measured 
from wea k tens ing surveys (|van Waerbecke et all (|200lf ): 
lRefregie"il (|2003|)1 and the Lyman - q fo rest absorption 
jCroft et all (| 19991 ): [McDonald et all (|2005l ')') are evidences 



toward an apparently clear picture of our universe at present. 
It is characterised by: I. spatial fiatness, II. a subcritical mat- 
ter content, HI. and accelerated expansion. 

But the clearness of this sketch poses a more interesting 
and deeper problem: how can we interpret all these features 
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models o f dark energy and dark ma tt er to alternative gravity 
theories (IPeebles fc Rathra I (120031): IPadmanabhan I (2003h : 
ICopeland. SamL and Tsuiikawal (200d )l" 

Unfortunately, so many data are not yet able to give us 
a definitive answer about the origin and nature of the ac- 
celeration of the Universe; not able to solve the coincidence 
problem, (namely, why we are currently detecting an energy- 
density for dark energy which is quite equal to the one of 
dust matter); and not able to state what the right EoS of 
dark energy is. 

Regarding the EoS, using current data one will be typ- 
ically only able to infer that the dark energy effective EoS 
parameter is close to —1. But any small deviation from 
this value could give a different theoretical scenario: if it is 
exactly equal (when of course referring to observational er- 
rors) to —1, we have a cosmological constant; if it is larger 
than —1, we have a quintessence model; while if it is smaller 
than —1 we have the so called phantom dark energy. In ad- 
dition, the data seem to indicate that the fractional energy 
densities of two main components of the Universe, i.e. dark 
matter and dark energy, are very similar at present, and the 
label "coincidence problem" has been coined to refer to this 
striking similarity. 

One of the most interesting solutions proposed to try 
and thr ow some l i ght on these questi o ns is o s cillati ng dark 
energy (|Dodelsonl ((20001 ): iFend ((2OO6I '): ISahnl (|2000l )). Such 
a model can easily solve the coincidence problem in a very 
natural way due to periods of acceleration, and can be also 
used as a good candidate for the unification of the late time 
acceleration (the one observed at present) with inflation (an 
early time acceleration period). In this context we have to 
underline the difference between assessing a periodic or non- 
monotonic potential and an oscillating dark energy EoS. In 
many cases such potentials do not give rise to a periodic 
w; one example can be found in iFriemanI (|l995h . where the 
proposed field is clearly periodic but the derived w can be 
well describe d by the GPL par amet rization for da rk energy 
introduced in IChevallieJ (|200ll ) and lLindeJ (|2003l ). 

In this p aper we are going to follow the method from 
iLindeij (|2006l '). by examining some directly proposed phe- 
nomenological periodic equations of state for the dark en- 
ergy using different cosmological observations. Specifically, 
we are going to set constraints on the location of the cen- 
tre of the range about which w oscillates and the amplitude 
of the oscillations, and we will also constrain the fractional 
energy density of matter. 

Actually, in the models to be considered there is an- 
other important parameter, the frequency of the oscillation. 
Relevant though it is, leaving this parameter completely free 
leads to a high dimensionality statistical problem, but given 
the precision of the data available at present it seems that 
problems of that sort cannot be tackled satisfactorily; i.e. 
there seems to be not enough quality in th e data t o con - 
strain more tha n two dark energy parameters iLinderl (|2005h : 
ISuUivarJ (|2008l ). In the literature on oscillating dark energy, 
the usual practice has been choosing a specific single value 
of this oscillation frequency and stick with it. I n orde r to fix 
it one may resort to an argument by iLindeJ (|2006l ) which 



^ This is the usual factor which relates pressure and density of 
any given component through the relation wx = Px / Px ■ 



suggests the lowest bound on this frequency for the data in 
use to be able to discriminate an oscillating behaviour from 
a monotonic one. We wish to carry out here a more thorough 
study of oscillating dark energy than previous works, and to 
this end we choose a discrete set of values of the frequency 
(above and below that bound) and then obtain constraints 
on the rest of parameters. Relaxing assumptions on the fre- 
quency, as compared to previous works, will allow us able to 
draw stronger conclusions; yet this is not the only novelty 
of our analysis, as will be shown in what follows. 

In this exploration of possible oscillating patterns in the 
expansion of the Universe induced by the dark energy com- 
ponent we also find it interesting to consider non-periodic 
oscillating models. Specifically we deal with scenarios which 
display a similar start off to the popular oscillating model 
by hinder, but then depart from it as their amplitude gets 
smaller as 2 grows and the distance between nodes tends 
to converge to a specific value. In this direction, we pro- 
pose a pair of models in which dark energy oscillations are 
modelled via special functions. Comparison between mod- 
els of that sort and the usual trigonometric parametrization 
provides hints about which features in the oscillations are 
favoured or disfavoured. 

Moreover, the present paper innovates in another di- 
rection: we choose a combination of datasets with inter- 
esting characteristics: we combine the statistically most 
powerful dataset available, the luminosity measurements of 
SNela, with other datasets: the luminosity measurements of 
Gamma Ray Bursts (GRBs) and direct Hubble data. GRBs 
are particularly useful for the study of oscillating models, as 
they typically inform us of higher redshifts than supernovae 
data, they improve the capability for detecting oscillating 
features (if they exist) at lower frequencies. In addition, the 
inclusion of the direct Hubble data can in principle enhance 
the sensitivity of our tests to the presence of oscillations, as 
the use of these data does not involve performing an averag- 
ing of the inverse of the Hubble factor, and then a possible 
smoothing of the oscillations is partially compensated for. 
Another point in favour of the usage of this particular com- 
bination of three datasets is the rather good concordance 
among them (see Fig. [!}, which applies at least for the case 
of a constant w quintessence, and therefore seems a priori 
a property that will be shared by models with a dynamical 
EoS parameter. 



2 OSCILLATING DARK ENERGY 
2.1 Our parametrizations 

We have just presented motivations for studying oscillat- 
ing dark energy. In general, in order to infer conclusions 
about the dynamical behaviour of dark energy and its con- 
sequences in the expansion of the Universe, one has to make 
some concessions to try and make the best out of a col- 
lection of noisy scattered data. One of the most popular is 
to propose a parametric reconstruction obeying some basic 
requirements. Our proposal fits precisely in this kind of ap- 
proach, and in order to examine the adequacy of oscillating 
patterns in the dark energy EoS we consider a simple and 
periodic phenomenologi cal parametriz ation for the EoS, as 
proposed and studied in iLindej (|2006l ) : 
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(a) ACDM (b) Constant EoS 

Figure 1. Credible intervals of the different observational data sets: green ones correspond to GRBs, purple to Hubble parameter data, 
orange to SNela and blue to the combination of all observational data. To construct these contour plots we have considered a prior on 
Cm and uiQ from WMAP-5years. 



w{a) = Wc — Asm{B\Tia — Q). (1) 
Equation ([1]) describes evolutionary dark energy with 

• © being the phase of the oscillation, which for simplicity 
is assumed to be zero. 

• Wc being the centre of the range over which w{a) oscil- 
lates. In the = case the parameter Wc is equal to the 
present value of the dark energy EoS, w{a = 1) = Wc = wq. 

• A being the amplitude of the oscillations, which obvi- 
ously must be non zero for a dynamical w. If this EoS is 
the effective realisation of a canonical scalar field, A should 
fulfil the constrain wq — A^ —1 ; but we leave it completely 
free and let the observational data speak out its value. 

• i? is the frequency of the oscillations; and there are 
some key remarks to be made about it. In order to notice 
distinctly the presence of an oscillatory behaviour B should 
fulfil the condition 



IB In a„ 



> 27r, 



(2) 



where amin is given by the highest redshift in the dataset. 
For our observational data this value is associated with a 
GRB at z — 5.6, for which we infer the constraint B > 3.3. 
Nevertheless, as SNela represent by far the dataset with 
the largest statistical power in our analysis and they span a 
smaller redshift range, it seems reasonable that an oscillating 
pattern will be only detectable for larger values of B. 

Previous works have attempted to set constraints in 
these parameters but they have not been able to distinguish 
an oscillatory model from one with a constant EoS. This is 
partly due to the data sets employed, but also because of 
excessive restrictions on the parameters. 

In this new attempt at exploring oscillating dark en- 
ergy we use a new combination of data sets so as to obtain 
more reliable constraints, but we also make a few other key 
changes: specifically we carry out an analysis for a discrete 
set of fixed B values to try and avoid the arbitrariness in the 
choice of this parameter present in works by other authors. 

In addition, and given the lack of grounds for very 
strong restrictions on phenomenological parametrizations of 
dark energy we have analysed two more models which rep- 



resent a slight departure from the perhaps excessively nicely 
shaped trigonometric parametrization by Linder. Our first 
proposal in this direction is 



A Ji(Blna-O), 



(3) 



where again we have set the phase O = so that Wc = wo, 
and Jn is the Bessel function of the first kind with n = 1. The 
Bessel function Ji{x) shows oscillations which are damped 
with growing x, as opposed to the constant amplitude of 
the trigonometric case. So in this case the EoS parameter 
w would have an oscillating trend modulated by a damping 
effect which makes the amplitude smaller and smaller as we 
go back in time (in the redshift space). 

As we have to set observational constraints on parame- 
ters, we have to compute the Hubble factor 



H'^{a) = f7ma ^ + (1 — ^m)a 



-3(l + «7{a)) 



(4) 



which depends on the present value of the fractional matter 
density fim, and on the averaged equation of state parameter 



■w{a) 



/,'"^°^^(a)dln(«) 
ln(a) 



(5) 



So, we will work at some points with this expression 
of the EoS parameter instead of the usual one. It can be 
noticed that for the two parametrizations presented above, 
and we realised that the averaged forms follow a common 
pattern; the parametrization given by the Eq. [1] takes the 
form 



■w{a) — Wo + 



A cos{B In a) - 1 



B ln(a) 
whereas the parametrization given by Eq|3] 
A Jo(Blna) - 1 



w{a) = Wo -I- 



B 



ln(a) 



(6) 



(7) 



With Eq. [7| and Eq. [6] as a reference, we propose another 
parametrization with a damped oscillating behaviour. 



10(a) =^10 + 



A {Tv/2)H-i{Blna) - 1 



B 



ln(a) 



(8) 
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where wq is the value of the EoS at present and -ff_i is the 
Struve function, Ha, of order a = — 1. The Struve function 
lies between the trigonometric case, with constant ampli- 
tude, and the Bessel function, being less acutely damped 
than the Bessel function. It can be checked that the 'w{a) 
function for this new proposal is an oscillating one as well, 
so this new model fits is the aim of the discussion. 

Anyway, we have to underline that our three w mod- 
els depend on In a (i.e. ln(l -|- z)) so these effects are quite 
smoothed and not so evident when depicted on even a large 
redshift interval such as the one limited by the re-ionisation 
one {z = 1089). Howe ver, such parametrizations, as first 
shown bv lLindeij (|2006l ). allow for analytical expressions of 
H{z) (see Eq. [5]) and thus result convenient for evaluation 
purposes. 

Interestingly, as we will discuss later, the main conse- 
quences of an oscillating equation of state are more evident 
on the deceleration q{z) function than on the Hubble func- 
tion or the dark energy fractional density ilx- 

2.2 Earlier Works 



In lRiesd (|2007l ) an interesting observational clue of a possible 
oscillating behaviour in the EoS was found. Fitting a quartic 
polynomial of w{z) to SNela observations indications that 
such oscillations might be really present were found. Even 
though the redshift range where these oscillations seem to 
be present is actually rather small; these results motivated 
a pleth ora of work s trying to analyse this possibility deeper. 
In 1x3 (|2005l 'l the formulation 



io(z) — Wo + wi sin{z), 



(9) 



was used to model a quintom scenario phenomenologically. 
The analysis was conducted using luminosity distances of 
SNela, the CMB shift parameter and the linear growth rate 
from large scale structure. Combining all these data sets 
the best fit values found were liio = —1.33, wi — 1.47. 
Such oscillating model differs very little from the linear one 
{w{z) = wo + wiz) when considering constraints from SNela 
only, and the oscillating case is only mildly preferred. This 
was somehow to be expected as for low redshifts the two 
parametrizations are very similar. Note that as in this study 
the oscillation frequency was fixed to a specific value we meet 

again the difficulties to make general infe renc es. 

An extension of the latter was done in I3cial (|2006l ). where 



10 (a) = Wo + wi sin (w2 In a + M3) , 



(10) 



was proposed, therefore now including also a period and a 
phase. One further degree of sophistication was introduced 
in the analysis as it was carried out using a modified version 
of CosmoMC so as to have into account dark energy per- 
turbations. The results indicate that wo is acceptably well 
constrained, while wi is poorly constrained, and W2 and W3 
can be regarded as completely unconstrained. This is not 
surprising, though, as highly dimensional parametrizations 
of dark energy seem to suffer this problem typically given 
the pre cision of ob servational data at present. 

In IXial (|2008l ) constraints on the cosmological param- 
eters of some oscillating models were analysed using simu- 
lated data from future Planck measurements. Two different 
parametrizations of the dark energy EoS were used, the CPL 
one 



w{a) = uiQ + Wi (1 — a) (11) 
and the a reformulation of the lLindeil (|2006l ) oscillating one. 



'j{a) — Wo + wi sin(TO2 In (a)), 



(12) 



having fixed the phase parameter to zero. In addition, one 
further choice was made: for simplicity and for focusing the 
analysis on low redshifts, the period parameter was fixed 
as W2 — 3n/2 (because it corresponds to a period in the 
redshift range < 2 < 2 where data from supernovae are 
more robust). Then they use present data from WMAP3, 
ESSENCE SNela, SDSS (for details see reference) to derive 
best fit values for the EoS parameters which are assumed 
as the fiducial values for deriving constraints from Planck 
mission details. The following results were obtained: 

, riq+0-S'i2+o.78i f-,n\ 

i.UO_o.587_i.57o Uoj 



wo 



-, ^0+0.15+0.36 
-i.UJ_o.l5_o.26i ^1 



for the CPL EoS, and 

n noi +0.320+0.534 
^0 = —0.9ol_o. 340-0. 7481 



Wl 



r, nco+0. 561+1. 037 
-U.Ob»„o. 591-1. 245 (,J^4j 



for the oscillating model. 

Assuming the CPL EoS they are able to conclude that 
a Quintom scenario with a. w{z) which crosses the phantom 
divide is favoured against the ACDM scenario (it is mildly 
favoured with present data while future surveys can provide 
narrower constrains) . Assuming the oscillating EoS the same 
conclusion can be derived, even if the constraints are weaker, 
and the low value for the amplitude of oscillation makes the 
constant w{z) (i.e. ACDM model) equally possible. 

For what co ncerns the oscillating model, these results 
are backed up in [lIuI l|2009l) , where the following 68% CL 
values were obtained: 



wo 



-0.958 ± 0.098, Wa 



0.030t°:r30 



(15) 

From these results one can conclude that the EoS fit with 
current data provides as good as ACDM, and that the os- 
cillation amplitude is limited: \wamp\ < 0.232 at 95%CL. 

An oscillating mod el dif f erent from the one proposed 
by Lind er is analyse d in ljaini (|2007l ) (inspired by a previous 
idea in iFend (|2006l ')). The authors used SNela, the CMB 
shift parameter and the measurement of the BAG peak from 
SDSS in order to constrain the sine version of the oscillating 
EoS written in the form 

w{a) — Wo — Asm{Blna) (16) 
and the alternative 

TO (a) = — cos (6 In a). (17) 

The period parameter b in the second model drives the ac- 
celerating/decelerating epochs, and in the limit of small b 
one has w{a) ~ —1. 

Final results show that b — 0.06 ± 0.01 at la level, 
a value which is consistent with limits required by CMB 
and correct power spectrum and which represents a clear 
evidence for an oscillating behaviour with a very long period 
(i.e. small but clear deviation fr om cosmologi cal constant). 

Another interesting work is iKurelj (|2008l ). where many 
different EoS were compared with a Bayesian approach. 
They worked with the so called linear EoS, i.e. CPL 
model; with pure oscillating models, the cosine one, w{z) = 
wocos(wcln(l + z)), and the sine one, w{z) = —1 -|- 
lOo sin(u;s ln(l + z)); with a damped version of previous os- 
cillating models, w{z) = wo{l + z)^ cos(uic ln(l + z)) or 
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■w{z) — —1 + wo(l + 2)^ sin('u;3 ln(l + z))\ and with a more 
complicated version of the dark energy EoS directly derived 
from the dynamics of the phantom scalar field (see reference 
for details). Analysing the values of the Bayesian evidences 
obtained fitting data (SNela, CMB, BAO) for each partic- 
ular model, ari d comparing them with the Jeffreys scale 
jjeffrevd (|l96ll )). they concluded that there is a substan- 
tial evidence for preferring pure oscillating model (sine) and 
the one derived from phantom field dynamics over the linear 
one, while damped versions are completely excluded. In ad- 
dition, they found no strong evidence to favour ACDM over 
the oscillating models. 

As stated in the Introduction, completely differ e nt ap - 
proaches are possible as well; for example, in ISaej (|2009l ') 
the starting point was an oscillating model for H{z) instead 
of wlz). Finally, it is possible to consider oscillati ng fields, 
which do not necessarily produce an oscillating EoS lFriemanI 
1 19951 ). 



3 OBSERVATIONAL DATA 

We have tested the possible periodicity of the Hubble func- 
tion by using three different observational data sets, i.e.: 

• the reconstructed Ifubble data given in ISternI (|2009l '): 

• the Con stitution Supernovae data set described in 
iHickenI (|2009D : 

• the Gamma Ray bursts lum inosity distances measured 
and analysed in iKodamal (|2008l ). 



3.1 Hubble parameter: Stern et al. 2009 data set 

Recently in ISterrJ (|2009l ) an update of the Ifubble func- 
tion H[z) data extracted from differenti al ages o f pass ively 
evolving galaxies previously published in lSimoi] (|2005l ) was 
presented. Constraining the background evolution of the 
Universe using these data is interesting for several reasons. 
Firstly, they can be used together with other cosmological 
tests in order to get useful consistency checks or tighter con- 
straints on models. Secondly, and more importantly, in con- 
trast to standard candle luminosity distances or standard 
ruler angular diameter distances, the Hubble function is not 
integrated over. This is a key point because if a periodic be- 
haviour is present in w{a) it should be directly detectable in 
H{a) while it could be lost in luminosity or angular diameter 
distances because of integration stages. 

The Hubble parameter dependence on the differential 
age of the Universe in terms of redshift is given by 

1 dz 



H{z) = - 



1 + z dt' 



(18) 



Thus, H{z) can be determ i ned fr o m measurements of dt/dz. 
As re porte d in [jimenej (j2002l l. Ijimenej lj2003l l. ISimonI 
l|2005h and [ster3l2Q0l)7 values of dt/dz can be computed 
using absolute ages of passively evolvin g galax i es. 

The galaxy spectral data used by SternI (|2009l ) come 
from observations of bright cluster galaxies done with 
the Keck/LRIS instrumenllj. The purposely planned Keck- 



survey observations have been extended with other datasets: 
SDSS improvements in calibration available in the Pub- 
lic Data Relea se 6 (DR6) have been applied to data in 
Ijimeriej (120031'): the SPICES infrared-selected galaxies sam- 
ple in lSternl (l2000l): and the VVDIS survey by the VLT/ESO 
telescope in lLe Fevrd (|2005l ). 

The authors of these references bin together galaxies 
with a redshift separation which is small enough so that the 
galaxies in the bin have roughly the same age; then, they 
calculate age differences between bins which have a small 
age difference which is at the sam e time larger than the 
error in the age itself (|Sternl \200^ ). The outcome of this 
process is a set of 11 values of the Hubble parameter versus 
redshift. A particularly nice feature of this test is that the 
sensitivity of differential ages to systemat i c erro r is lower 
than in the case of absolute ages (|jimenezl (|2004l )). 

Observed values of H{z) can be used to estimate DE 
parameters by minimising the quantity 



{Hiz,-{e^)}-H,U^j)f 



(19) 



where Hq = 100 h wiU be fixed as h = 0.742 (|Riesj 
(|2009l )). while the vector of model parameters, 6i, will be 
Oi = (ft Wo, A, B). 



3.2 Supernovae: Hicken et al. 2009 data set 

We use one of the most re cent SN e la sam ples, the Consti- 
tution sample described in iHickenI (|2009l ). which is a data 
set obtained by combining the Union data set by iKowalskil 
(2008) with new 90 nearby objects from the CfA3 release 
described in IHickenI (j2009Al ). 

The Union SNe compilation is a data set of low-redshift 
nearby-Hubble-flow SNe and is built with new analysis pro- 
cedures for working with several heterogeneous SNela com- 
pilations. It includes 13 independent sets with SNe from the 
SCP, High-z Supernovae Search (HZSNS) team, Supernovae 
Legacy Survey and ESSENCE Survey, the older data sets, as 
well as the recently extended data set of distant supernovae 
observed with HST. After various selection cuts were applied 
in order to create a homogeneous and high-signal-to-noise 
data set, a final collection of 307 SNela events distributed 
over the redshift interval 0.15 z 1.55 was obtained. 

The CfA3 data set was originally made of 185 multi- 
band optical SNela light curves obtained at the F.L. Whip- 
ple Observatory of the Harvard-Smithsonian Center for As- 
trophysics (C fA); 90 of the or iginal 185 objects passed the 
quality cuts of lKowalskil (|2008l ) and were added to the Union 
data set to form the Constitution one. 

The statistical analysis of the Constitution SNe sample 
rests on the definition of the modulus distance, 

m(2,) = 51ogio[di(z,, {6(,})] + fio, (20) 

where dL{zj, {6i}) is the Hubble free luminosity distance 

_1 



dLiz,{e,}) = (l + z) / dz'—j 



(21) 



bee ISterni ||2009B|) for a detailed description of the observations, 
reductions and the catalog of all the measured redshifts 



The best fits to be presented will be obtained by minimising 
the quantity 
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where the are the measurement variances. The nuisance 
parameter fio encodes the Hubble parameter and the abso- 
lute magnitude M, and has to be marginalised over. Giv- 
ing the heterogeneous origin o f the Con s tituti on d ata set, 
and th e procedures described in iKowalskil (|2008h and lHickenI 
l|2009l ) for reducing data, we have worked with an alterna- 
tive version of Eq. (|22|) . which consists in minimizing the 
quantity 

2 ci 
XSN{{di}) = ci - — 
C3 

with respect to the other parameters. Here 



(23) 



307 

307 

3 = 1 



{^{zj;^io = 0, {Oi)} - fiobs{zj)f 



{n{zj;fio = 0, {Oi)} - fiobsizj)) 



307 



C3 



(24) 



(25) 



(26) 



It is trivial to see that Xsjv is just a version of XsiV) min- 
imised with respect to To that end it suffices to notice 
that 



XSn(MO. {^i}) = Cl - 2C2^iO + C3/io 



(27) 



which clearly becomes minimum for /lo = C2/C3, and so we 
can see xsn = xiN(Mo ~ 0, {Si}). Furthermore, one can 
check that the difference between Xsiv s-^d Xsiv is negligible. 



3.3 GRBs: Kodama et al. 2008 data set 



The GRBs sample, described in lKodamal (|2008l ). is made of 
33 GRBs within the redshift interval z < 1.62 and 30 GRBs 
in the redshift interval 1.8 < z < 5.6. It is well known that 
GRBs are not standard candles as SNela; at the same time 
they contain a lot of important information about high red- 
shift properties of the universe which cannot be derived from 
SNela. So their combined use can bring important and com- 
plementary information about the reconstruction of dark en- 
ergy and gives us the possibility to detect eventually traces 
of an oscillatory behaviour on a larger redshift range. The 
calibration of GRBs data can be done in several ways, and 
many empirical formulas have been given f or describ i ng the 
peak energy-peak luminosity correlation. In lKodamal (|2008l ) 
the peak energy-peak luminosity correlation is described by 
the so called Yonetoku relatioi^ of the GRBs sub-sample in 
redshift interval z < 1.62 is calibrated without assuming any 
cosmological model and using the luminosity distance of the 
objects considered it is estimated from those of SNela with 
redshift z < 1.755. The calibrated Yonetoku relation is then 
finally applied to the high redshift GRBs sub-sample with 
redshift in the interval 1.8 < z < 5.6. Final data consist of a 
set of calibrated luminosity distances, so that we can define 
the contribution to the total chi-square as: 

3 See lYonetokul b004] 



^2 (fni^ ^{dL{zr,m}-driz,)f 



3 = 1 



< (-3) 



(28) 



where the cr| are the measurement variances. 



4 STATISTICS AND DATA ANALYSIS 

We will explore the probability distributions of our problem 
with Markov Chain Monte C arlo (MCMC ) methods . Thes e 
methods (f ully described in iBerd (|2004l '). iMacKavl (|2002h . 
iNeall (|l993l ) and references therein) extract samples sequen- 
tially using a probabilistic algorithm (we have chosen the 
Metropolis-Hastings algorithm which we describe below) 
based on the Bayesian statistical approach. The main 
problem when running a Markov chain concerns how to 
avoid biased inferences or underestimations on errors of 
the theoretical parameters. The problem of the convergence 
is the main one among the problems of that sort, and it 
can be formulated as follows: how can one be sure that 
the properties of a sample from the MCMC algorithm are 
a good representation of the unknown distribution to be 
explored? In the following we discuss the details algorithm 
and the solution adopted for convergence. 



4.1 Markov Chain algorithm and convergence test 

The Metropolis-Hastings algorithm works as follows: start- 
ing from an initial parameter vector p (in our case p = 
{^IjnjWQ, A, B)), one generates a new trial point p' from a 
proposal density q{p',p), which represents the conditional 
probability to have p' given p. This new point is accepted 
with probability 



, ,, . ( P p' d g p',p 
nP> P = mm <^ 1, ) / ^ 



(29) 



where P(p|d) is the conditional probability to have the pa- 
rameter set p given the observational data d. Then, from 
Bayes' theorem it follows that this probability is: 



P(p|d) cc i(dlp)P(p) 



(30) 



with L(d|p') being the likelihood function containing infor- 
mation from the data, and P(p) the prior on the parameters 
which contains all the supposed information on them before 
observing the data. If the chain moves to the new set p', 
then one says that it has been accepted, otherwise it has 
been rejected. 

We choose the prior depending on the physical require- 
ments a given parameter has. The proposal density is typ- 
ically a Gaussian symmetric with respect to the two vec- 
tors p and p', namely g(p, p') oc exp(— Ap^/2(Ty ), with 
Ap — p — p', so that, from the detailed balance equation, 
we know that the final probability distribution is stationary 
under the Markov process. An important quantity in the 
testing proposal distribution is its dispersion ctt; as we will 
discuss below, we decide not to take a fixed value for it, but 
instead, we let it depend on the value of the parameters at 
any step. 

One can say that a chain has reached convergence when 
the statistical properties of the extracted samples can de- 
scribe the statistical properties of the unknown probability 
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distribution with good accuracy. Probability theory says that 
Markov processes will reach the exact final distribution in an 
asymptotic way, which means a sample's length will become 
infinite in an infinite computation time. Of course, as one 
is forced to operate with finite length samples, the question 
arises of how this truncation can bias the final statistical 
results, and if there are any parameters that can be able 
to give information about the goodness of the process. In 
the literature the mos t used paramet er for this task is the 
convergence ratio (see lDunklevI (|2005l 'l'). defined as 



(31) 



This is the ratio between the variance of the mean of the 
samples and the variance of the underlying unknown dis- 
tribution (we will operate with standard distributions so 
that (Tg = 1). Then r is required to be below a cut-off 
limit value, typically 0.01, to have a guaranteed convergence 
of the chain. This parameter is used in other convergence 
tests, such as the Gelmann-Rubin test (see lGelman fc RubinI 
which runs many parallel multiple chains and esti- 
mates r at any step; but we regard it a time and hardware 
expensive test. 

An alternative solution to this probl e m in the spectral 
analysis approach proposed bv iDunklevI (|2005l '). Is it clear 
that all the steps in MCMC are correlated; this correlation 
is somewhat intrinsic to the code, at least before having 
reached convergence, but it depends also on the value of 
(Tilfl. But what is even more important is the behaviour of 
the correlation when the convergence has been reached: in 
this case the MCMC will sample from the underlying distri- 
bution and it will work like a random sampler, so that there 
will be no correlations in this regime. 

This is the key idea of the test bv iDunklevI (|2005l '): if 
we take the power spectra of the MCMC samples, we will 
have a large correlation on small scales, but the spectrum 
will become flat (like a white noise spectrum) when conver- 
gence has been reached. Then, checking the spectrum of just 
one chain (instead of many parallel chains as in Gelmann- 
Rubin's test) will be sufficient to assess that convergence 
has indeed been reached. We will give just a short account 
of the steps to be followed to implement t he test, but for a 
more detailed reference see IDunklevI (|2005l '). 

In brief, we calculate the discrete power spectrum of the 
chains. 



or in equivalent logarithmic form: 



with 



1 ""'^ 



exp 



.2-Kj 



(32) 



(33) 



where A'' and are respectively the length and a given 
element of the sample from the MCMC, j = 1, . . . , N/2 - 1, 
and the wave number kj of the spectrum is related to the 



index j by the relation 
analytical template: 



P{k) 



'i + {k*/kY 



211 j /N . Then we fit it to an 



(34) 



In P, = In Pa + In 



(35) 



where 7 = 0.57216 is the Euler-Mascheroni number and 
rj are random measurement errors with < rj and 
< riTj >= Sijit^ /&. The fit provides estimates of three pa- 
rameters, but only two of them are fundamental to our anal- 
ysis. The first one is Pq, which is the value of the power 
spectrum extrapolated for fc ^ 0; this is an important pa- 
rameter because from it we can derive the convergence ratio 
using r ~ Po/N . The second important parameter is j* (the 
index corresponding to fe*), which is related to the turnover 
point from a power to a flat spectrum; the estimated value 
of j* has to be >> 20, so one can be sure that the num- 
ber of points in the sample coming from the convergence 
region are larger than the number of noisy points. If these 
two conditions are met for all the parameters, then the chain 
has reached convergence, and the statistics from the MCMC 
procedure describes well the u nderlying probability distri- 
bution. Following the advise in iDunklev (|2005h we perform 



the fit over the range 1 ^ j ^ j„ 



with 



lOj*, 



where a first estimation of j* can be obtained from a fit 
with jmax = 1000, and then perform a second (or even a 
third) iteration to have a better estimation of it. 



4.2 Model selection tests 

After having estimated the value of the parameters set by 
using the MCMC approach, we need a tool for comparing, 
selecting and testing the statistical goodness of our results. 
The related literature is too extensive for being reviewed 
here, so we will only refer in some detail to the tools we 
have chosen. 

Since the MCMC technique is based on a Bayesian ap- 
proach, the best way for comparing models is arguably the 
Bayesian Evidence. It is defined as 



E = / c{e)P{e)Ae, 



(36) 



where C is the likelihood function, 6 is the parameters vector 
and P{9) is the prior distribution for the parameters. It is 
clear from its definition, that the evidence of a model is the 
average likelihood of the model with respect to the prior: 
models which flt the data well and make narrow predictions 
are likely to fit well over much of their available parameter 
space, giving a high evidence. So using evidence for com- 
paring models is a very appropriate task. Model comparison 
requires defining then Bayes factor. 



E{Mi) 
E{M,) 



(37) 



^ For a detailed analysis of this aspect see IDunklevI ||2005| ) 



which is the ratio between the evidence values of two models. 
Mi and Mj. If Bij > 1 then the model Mi is preferred 
with respect to the model Mj. B y conven ti on, B ayes factor 
is judged on the Jeffreys' scale (|jeffrevsl (|l96lh l: for 1 < 
In Bij < 2.5 there is a "substantial" evidence in favour of the 
model with the greatest Bayesian evidence; for In Bij > 5 
the evidence is "decisive" . 

Moreover, we have to underline in favour of the 
Bayesian evidence, that it is a full implementation of 
Bayesian inference and can be directly interpreted in terms 
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of model probabilities. Unfortunately, being a highly-peaked 
multi-dimensional integral, its estimation typically requires 
a hard and challenging numerical effort. Even if some algo- 
rithms have been found which simplify this operation, it is 
always preferable to have easier tools for estimating it and 
comparing models. 

The easiest and most used tools are different versions 
of the Information Criteria. Generally, the introduction of a 
higher number of parameters improves the fit to the chosen 
dataset, regardless of whether or not these new parameters 
are really relevant. As a consequence, the simple compar- 
ison of the maximum likelihood value of different models 
will tend to favour the model with the highest number of 
parameters. The information criteria work just in this direc- 
tion: they compensate this behaviour by penalising models 
which have more parameters. 

The first test is the Akaike Information Criterion (AIC) 
defined as 



AIC 



-2\nC + 2k 



(38) 



where C is the maximum likelihood value a nd fc is the 
number of parameters of the model (|Akaikel (|l974l )V The 
AIC is derived by an approximation of the KuUback- 
Leibler information entropy, which measures the difference 
between the true data distr ibutio n and th e mod el one 
IIBurnham fc Andersoril (|2002l ) and iTakeuchil (|2000l ll. The 
best model is the one which minimises the AIC, and no re- 
quirement for the models is asked for. 

There also exists an AI C version for small samp l e size s, 
the corrected AIC, AICc (|Burnham fc AndersonI (|2002l ')') 
given by 



AICc = AIC + 



2k{k + l) 



N ■ 



1' 



(39) 



where A'' is the number of points in the dataset. Since the 
correction term disappears for large sample sizes, A'^ >> k, 
we will use th is last definit ion for comparing models (as 
pointed out in iLiddlj (|2007l ). it is always preferable to use 
the corrected version rather than the original one). 

From the same principles of AIC (minimisation of the 
KuUback-Leibler information entropy) another comparing 
tool can be derived, the Residual Informat ion C r iterio n 
(RIC). We will use the corrected version in iLend (|2007t ). 
where the RICc is defined as 



RICc = -21n£ + (fc-l) + 



4k 



N-k- 



(40) 



When N » k, RICc has a smaller penalty than AICc. 

Generally, when A AIC > 1, it follows that the two mod- 
els are significantly different, and the one with the lowest 
value of AlCid the preferred one. Finally, we have to keep in 
mind that AIC tends to favour models with a high number of 
parameters and it is "dimensionally inconsistent", namely, 
that even as the dataset size tends to infinity, the proba- 
bility of the AIC incorrectly picking an ov erparametrized 
model does not tend to zero (|Liddlel (|2004 ) and references 
therein). Let us thus consider alternatives. 

The Bayesian Infor mation criterion (BIC) was intro- 
duced in ISchwarj (| 19781 ') and is defined as 



BIC 



-2ln£ + klnN 



(41) 



Again in this case the best model has the lowest BIC, and it 
is clear from this expression that BIC penalises models with 



a high number of parameters more than AIC. Being the BIC 
a good approximation for twice the log of the Bayes factor, 
it can be compared with Jeffreys' scale, so that ABIC > 5 
means a "strong" evidence in favour of the model with lowest 
BIC values, while for ABIC > 10 this evidence is "decisive" . 

We will also use t he D eviance information criterion 
(DIG) of lSpiegelhaiteil (|2002D . which mixes elements from 
both Bayesian and information theory. It is well suited to 
our case because it is easily computable from posterior sam- 
ples such as those coming from MCMC runs. It relies on the 
definition of the effective number of parameters, pn, also 
known as the Bayesian complexity. It is defined as 



PD = D{9)~D{e), 
where 

D{e) = -21n£(6i) -fC 



(42) 



(43) 



with C a constant which vanishes from any derived quantity, 
and the chi-square defined as usual as = —2 In/!. This 
definition shows that pn can be considered as the mean de- 
viance minus the deviance of the means, and this is the key 
quantity in estimating the degrees of freedom of a test. Fi- 
nally the DIG is defined by 



DIG = D(e) + 2pD = D{e) + pD 



(44) 



where one can recognise a similar-to-AIC formulation in the 
first expression, while a Bayesian definition and measure of 
model adequacy (penalised by an additional term related 
to the model dimensionality) is implicit in the second one. 
The DIG is also useful for another reason: it overcomes the 
difficulties AIC and BIC have to discount parameters which 
are unconstrained by data (BIC is of even more suspicious 
validity when there is any parameter degeneracy). Finally, 
DIG (like BIG) is not dimensionally inconsistent, so it is able 
to detect wrong high dimensionality parametrizations. 



5 CONSTRAINTS AND ASSUMPTIONS 

We have implemented a few priors for running our MCMCs. 
The main one has been to set the control < f2m < 1, 
which is a minimum physicality requirement. We have also 
set mild Gaussian priors on fi™, and wQ with the 3cr error 
bar from WMAP5 as a reference. 

As we stated in the previous sections it would be possi- 
ble to set physical limits on the parameters of the oscillating 
model, such as A and B. The amplitude A should have a 
value which depends on the theoretical scenario chosen to 
be the fiducial one. In the case of ACDM, one could require 
that the minimum value for the EoS was too — A ^ —1. But 
it is clear that with this case will exclude the possibility of 
a phantom behaviour, which is a possibility present data do 
not rule at and in some cases seem to be the preferred one. 
For that reason regard leaving A free as the best option. 

At the same time, the frequency (period) parameter B 
should be suject to | _B In ttmin| > 2-k. The highest redshift 
of our observational data corresponds to a GRB observation 
at z — 5.6. For that choice, B should be B > 3.3. Anyway, 
we have no a priori strong clues about the validity of a peri- 
odic oscillating EoS, so there is scarce guidance regarding a 
suitable lower bound on B. We could have oscillations with 
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a period bound given by the highest redshift in the super- 
novae data, so in this case we would have B > 5.7; or we 
could have no detectable oscillations at all in our redshift 
range so that a really small value of B could turn out to be 
the best fit. In addition, rigourosly the bound on B we are 
commenting about (proposed by Linder) makes only sense 
for the sine oscillating model, which has a valid definition 
for the oscillation period and not for other two models. Any- 
way, also in this cases the B parameters can be related to 
periodic properties of the dark EoS parameter, so we treated 
it on the same footing in the analysis of the three models. 

There is also another problem concerning the number 
of free parameters one introduces in a model. In our case 
we would have a dark energy EoS with three free param- 
eters, namely wq, A and B (which become four parame- 
ters because of the presence of Qrn in the expression for the 
Hubble function and luminosity distance) . This poses a well 
known problem in reconstructing or modeling the EoS with 
parametric relations: how many parameters can we enquire 
about? Following the Occar's razor prescription one could 
be tempted to choose the minimalistic option: models with 
few parameters are the best ones. But sometimes this could 
be a not physically good choice: complex systems could re- 
quire complex analytical formulas and a large number of pa- 
rameters could describe a more suitable behaviour of dark 
energy. It is also possible that not all the introduced param- 
eters are really free, and there is a correlation/dependence 
between some of them; but this cannot be known a priori 
when proposing a new model. 

The only solution is to decide depending on the phys- 
ical problem one has to face with. In preliminary runs we 
left all the parameters free, but it soon emerged that there 
is a strong degeneracy between some of them, mainly be- 
tween A and B. While Qm and wo were well constrained, 
the two main parameters of our oscillating models showed a 
degeneracy which made them eventually unconstrained and 
yielded no satisfactory information about our proposed EoS. 

So we turned to another way to proceed: we fixed the 
value of -B to a set of discrete values scanning entirely the 
range of values which could potentially lead to detectable 
oscillations of our observable functions given the criteria dis- 
cussed above. 
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Figure 2. In these figures we plot the redshift- variation of the 
EoS and the aeceleration for the best values of the parametriza- 
tions. The colors have the same meaning of Figs. (2). 
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6 RESULTS 

Tables [TJ [2] and [3] summarize our results. As ours is a more 



Table 2. Averaged EoS: w{a) = wq + 4 '^°'fn'("')'" 
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in Fig. [6] A fit of yl as a function of B using a linear relation 



exhaustive analysis than previous ones in the literature it is ^ 1.0806 53.39 0.25l°, p2 -O.lll^g™ -0.891°, " 

clear that we can draw stronger conclusions. The main one ^ 1.0801 19.09 0.25lo p2 -0.02to^l -0.9llo^3 

is that the current data do not seem to give as strong con- ^ 1.0801 10.02 0.26lo j^2 -O.OOlJ^JIJ -0.90to\\l 

straints on A as on the other parameters Qm and wo, but a 1.0818 6.39 0.261q p2 "O ^olo gg — 0.86lQ j^4 

clear trend in j4 can be guessed, which becomes quite evident "e" 1 0809 4.55 0.25_(, p2 ~019-o!55 ~0-8^-o.i7 
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turns out to be the best one for the sine model; whereas, 2 '^■^'^^'^ ^ '^ tS S? zS gi lS ?n 

f 13 1 J c+ f +• ..q +• 1 +• • ■¥ 1-0808 2.32 0.261°-''^ -0.58i°'^^ -0.7llg-i2 

for the Bessel and Struve function a quadratic relation is .,.+0:02 „ 7^+0 el „ k^+o- w 

preferred — 1.0807 1.98 0.25_o 02 -0-75_o^6 -U.65_o 23 

If we pay attention to the behaviour of wq, which is 
present value of EoS, we see that is very well constrained. 
There is a slight difference between the sine and the Bessel 
or Struve cases, but in all cases we can exclude phantom 
like behaviour at la level (the prior does not in principle 
hinder it as it rather weak). In Fig[6] we can see that the 
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Table 3. Averaged EoS: w{a) ■ 
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Figure 3. In these figures we plot the free parameters 
and the reduced value versus the fixed B value for the different 
models. The blue points refer to the sine model (Eq.(6)); the 
red ones to the Bessel model (Eq.(7)), and the green ones to 
the Struve model (Eq.(8)). The gray line in Fig. [6] and Fig. |6] 
corresponds to ACDM values. In Fig.[6]the gray region indicates 
the values of f!m within \a error that describe ACDM. 



behaviour for too as a function of B agrees with the tendency 
of the amphtude A: aa B grows the value of A becomes more 
negative so luo moves to less negative values, see EqO and 
Eq[T] If we explore by means of fits how w and B are related, 
we find the same pattern as for A, the linear fit is preferred 
for the EoS with a sine form, and the quadratic one for the 
EoS with the Bessel or Struve function. 

The remaining parameter, fim, is very well constrained 
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and fully agrees with literature expected values, being fim — 

0. 25 and this value changes negligibly with B, see Fig|6] 

So, with the current data it seems that B cannot be 
really constrained, that is, all B values seem to be of similar 
statistical validity. We have a very slight preference for val- 
ues which are different from those chosen by other authors. 
Generally the chosen value is B = 37r/2 as corresponds to 
the typical supernovae redshift range. In our analysis the 
lowest value of for chi-square does not correspond to that 
value of B\ for the sine and Struve models the minimum 
is for B = IItt/G which means that we can detect oscilla- 
tions within a redshift z ~ 1.98, comprised by our data. For 
the Bessel oscillating model we have a chi-square minimum 
value at B = 27r/3 which will need a redshift z ~ 19.09, 
which is outside of our observational data range. Fig. [2] re- 
fiects clearly these behaviours: the best values, from the sta- 
tistical point of view, for the sine and Struve models depend 
allow us to detect a complete oscillation in the range of our 
observational data. 

Nevertheless, focusing our attention on Fig. [6] we can 
see that for almost all the values of the frequency, B, the 
values of chi-square of the sine parametrization are smaller 
than those of others parametrizations. If we take this into 
account and the fact that the best value of this parametriza- 
tion corresponds to a i? = lln/6, we could say that obser- 
vational data shows a preference for a periodic EoS with a 
small period. 

At the same time, if we look at the variation of the 
acceleration parameter (Fig[6l) in the recent past we can 
observe that it has recently pea ked and is slow i ng do wn at 
present as have been pointed in IShafieloo et all (|2009l '). 

Moving to the statistical side of the analysis, we have 
to argue if the proposed models are reasonably good or not, 
and above all if they can compete with the the concordance 
ACDM model. From our results we conclude an oscillating 
pattern in the dark energy is an admissible possibility, as 
there is so far no concluding evidence against it; and from 
some statistical perspectives they even represent a better op- 
tion that its main competitor, ACDM. Tables [4] summarizes 
our findings on the statistical side. 

If we take a look at the reduced chi-square values, we 
can see that all the models have lower values than the ACDM 
one, even if the differences are really small. In particular, the 
sine model seems the best one. If we take a look to the AICc 
values, we see that the sine model is just in a border line 
position if we consider the limit we discussed in Section [4.21 

1. e. AAIC > 1. On the contrary, we should absolutely reject 
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Table 4. Statistical criteria 



Model 


x= 




AAICc 


ABIC 


ARICc 


ADICc 


ACDM 


507.7003 


1.0802 














Sinus 


504.6550 


1.0783 


0.9975 


9.2687 


-1.0282 


1.9145 


Bessel 


505.4669 


1.0801 


1.8094 


10.0806 


-0.2163 


2.0344 


Struve 


505.0848 


1.0792 


1.4273 


9.6985 


-0.5984 


2.0023 



the Bessel and Struve models. But the situation changes 
when considering RICc: we know that it has a small penalty 
than AICc when A'' >> fc as it is in our case, where we have 
A'^ = 472 and fc = 3. And we see that its values favour the 
oscillating models with respect of ACDM; they are even neg- 
ative, meaning that the RICc of oscillating parametrizations 
are better. 

The situation reverses again when moving to BIC; if 
we compare the values obtained with the Jeffreys' scale, we 
should conclude that there is an almost decisive evidence 
against oscillating patterns. But we have to remember that 
BIC has some problem when facing degeneracies between 
parameters, like the behaviours of the amplitude and of wq 
seem to reflect. 

These degeneracies would well be the reason of the "bad 
results" offered by BIC, which is also challanged by another 
criterion, DIC, which looks more favourable with oscillating 
dark energy. Since DIC relies also on the Bayesian approach, 
we can apply Jeffreys' scale to it in the same fashion as 
above, and from this we conclude there is not a significant 
evidence against oscillations. 



7 CONCLUSIONS 

In this paper we have performed a quite detailed analysis to- 
wards the detection of oscillating patterns in the dark energy 
EoS. We have considered different phenomenological models, 
starting from the original sine models and then introducing 
two new ones, based on the special functions. Those new 
models differ from the sine one mainly because they show a 
damped amplitude in the oscillations when moving into the 
past. 

Compared to prior works devoted to oscillating dark 
energy we can highlight that fact that instead of fixing the 
frequency parameter to a particular value, we have explored 
a discrete set of frequency values. 

We have also introduced novelty with respect of the 
datasets used, as we take direct measurements of the Hubble 
factor and GRB luminosity data, which are arguably two 
directions of improvement as the sensitivity to oscillating 
patterns gets increased and a wider redshift range comes into 
play. Our theoretical setup has been complemented with a 
detailed statistical analysis using different model selection 
tools. 

Numerical results show that while parameters like the 
matter content, Q.m and the present value of dark energy 
EoS, wq can be constrained very well, this is not true for 
the amplitude and the frequency. In particular, between the 
amplitude A and wo seems to be working a degeneracy that 
cannot be solved. About the frequency, we can say that chi- 
square values are not really in favour of a particular value, 
being all in a very narrow range. But the best values favour 



values of the frequency which mean oscillation detectable 
inside the present redshift range of SNela. 

The statistical analysis seems does not seem to provide 
such a conclusive answer as desirable, though we think that 
oscillations can be considered as a possible alternative to 
a ACDM model for addressing the well know problems it 
suffers from (as review in the Introduction). All but one 
(BIC) of the statistical criteria used considered here are keen 
to the possible detection of oscillations in the EoS, there is 
even one of them (RIC) which seems to favor an oscillating 
pattern against a cosmological constant. 
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